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Three different cluster expansions for the evaluation of correlated one-body properties of s-p and 
s-d shell nuclei are compared. Harmonic oscillator wave functions and Jastrow type correlations are 
used, while analytical expressions are obtained for the charge form factor, density distribution, and 
momentum distribution by truncating the expansions and using a standard Jastrow correlation func- 
tion /. The harmonic oscillator parameter b and the correlation parameter j3 have been determined 
by a least-squares fit to the experimental charge form factors in each case. The information entropy 
of nuclei in position-space (5V) and momentum-space (Sk) according to the three methods are also 
calculated. It is found that the larger the entropy sum S — S r + Sk (the information content of the 
system) the smaller the values of \ 2 - This indicates that S is a criterion of the quality of a given 
nuclear model, according to the maximum entropy principle. Only two exceptions to this rule, out 
of many cases examined, were found. Finally an analytic expression for the so-called "healing" or 
"wound" integrals is derived with the function / considered, for any state of the relative two-nucleon 
motion and their values in certain cases are computed and compared. 
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I. INTRODUCTION 



The effect of short-range correlations (SRC) to the one-body properties of nuclei is an old but challenging and 
appealing problem. In general, the account of SRC is important for the description of the mean values of some two- 
body operators, such as the ground state energy of nuclei but it is also of interest to investigate the SRC contribution 
to simpler nuclear quantities related to one-body operators such as the form factor (FF), density distribution (DD) 
and momentum distribution (MD). It has been shown that mean- field theories can not describe correctly MD and DD 
simultaneously M and the main features of MD depend little on the effective mean- field considered |2j . The reason is 
that MD is sensitive to short-range and tensor nucleon-nucleon correlations which are not included in the mean-field 
theories. We note however that the choice of a single particle potential having a short range repulsion could play a 
role in improving somehow the values of MD || . 

The experimental evidence obtained from inclusive and exclusive electron scattering on nuclei established the 
existence of a high- momentum component for momenta k > 2 fm _1 [Q-Q. It is well known, that the independent- 
particle model (IPM) fails to reproduce the high momentum transfer data from electron scattering in nuclei. That 
is, the IPM is inadequate to reproduce satisfactorily the diffraction minima of the charge FF for high values of 
momentum transfer. Therefore, although single-particle potentials of the type mentioned above, that is with a short- 
range repulsion lead to certain improvement, theoretical approaches which take into account SRC due to the character 
of the nucleon-nucleon forces at small distances, are necessary to be developed. 

In this effort, two main problems appear. The first one is the type of SRC which must be incorporated to the 
mean-field nucleon wave function and the second one is the type of cluster-expansion to be used which is connected 
with the number of simultaneously correlated nucleons. 

In the present work we consider central correlations of Jastrow type |8| while three different cluster expansions are 
considered. The first two types of expansions, named FIY (Factor, Iwamoto and Yamada) and FAHT (Factor, 
Aviles, Hartogh and Tolhoek) |nj respectively, were developed by Clark and co-workers lO,£2) while the third one 
named LOA (Low Order Approximation) was derived by Gaudin, Gillespie and Ripka fl3| , |14|| . 

The FIY expansion, truncated at the two-body terms, was used for the calculation of the charge FF and DD |u| 
and MD fuj in s-p and s-d shell nuclei while the LOA, truncated at the two-body terms and including a part of the 
three-body term, was used for the calculation of the above one-body quantities in the closed shell nuclei 4 He, 16 and 
40 Ca [17|] as well as of the bound-states overlap functions, separation energies and spectroscopic factors in 16 and 
40 Ca [lqj. The FAHT expansion, truncated at the two-body terms was used for the evaluation of the charge FF JToj ] 
and nuclear ground state energy of 4 He and 16 O [po| . In the present paper the FAHT expansion is used in addition 
for the evaluation of the FF, DD and MD in s-p and s-d shell nuclei. 
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The present work is, in a way, a generalization of Ref. |21| where a comparison of various cluster expansions for the 
calculation of the charge FF of 4 He was made. In this generalization, the above mentioned three types of expansions 
are applied and compared for the one-body characteristics of s-p and s-d shell nuclei. 

The comparison of the three truncated expansions can be made, as usually, by comparing \ 2 ( m computing the 
FF) i.e. the smaller the x 2 , the better the quality of the corresponding expansion. In the present work we introduce 
also an information-theoretical criterion in addition to \ 2 ■ Information-theoretical methods [^2|-[32"| play an important 
role for the study of quantum-many body systems. It has been found in Ref. ]3l| ] that interesting properties of the 
information entropy S hold for various systems. For instance, it was shown that S = a + blnN where N is the number 
of particles in nuclei, atomic clusters and atoms. In a previous work |30| it was found that the larger the S, the better 
the quality of the nuclear model. Here we apply this idea to compare various cluster expansions. It turns out that 
this is the case i.e. the larger the S the smaller the values of x 2 , for various nuclei and expansions, with only two 
exceptions. 

The paper is organized as follows. In Sec. II the general expressions of the one-body density matrix (OBDM) for 
the three types of expansions are given. Numerical results are reported and discussed in Sec. Ill, while the summary 
of the present work is given in Sec. IV. Finally, some details of the FAHT expansion as well as for the calculation of 
healing integrals are given in Appendix I and II, respectively. 



II. CORRELATED ONE-BODY PROPERTIES 



A. General definitions 



The key of the description of the one-body properties of nuclei is the OBDM p(r, r'), which for a system of A 
identical particles is defined J3^,Q in terms of the complete wave function V?(ri, T2, r^) by 

p{r,r') = y**(r,r 2 ,...,r A )*(r',r 2 ,...,r A )dr 2 ...dr J 4, (1) 

where the integration is carried out over the radius vectors and summation over spin and isospin variables is implied. 

In the case where the nuclear wave function ^(ri, r 2 , . . . , r^) can be expressed as a single Slater determinant 
depending on the single-particle wave functions we have 



WD (r,r')=X>:(r)&(r / ). (2) 

*=i 

The diagonal elements of the OBDM give the DD, p(r, r) = p(r) while the FF is the Fourier transform of it 

•F(q) = J exp[iqr]p(r)dr, (3) 

and the MD is given by a particular Fourier transform of the OBDM 

n(k) = ^ J exp[zk(r - r')]p(r, v')dvdv' . (4) 

The second moment of the DD is the mean square radius of the nucleus while the second moment of the MD is 
related to the mean kinetic energy. 

We also define the information entropy sum 

S = S r + S k , (5) 

where 



S r = - I p{r)\np{T)dv (6) 
is the information entropy in position-space and 

S k = - I n(k)mn(k)dk (7) 



is the information entropy in momentum-space. 

S is a measure of quantum-mechanical uncertainty and represents the information content of a probability distri- 
bution, in our case of the nuclear density and momentum distributions. In the present work, we employ in calculating 
S a normalization to the number of particles A for p(r) and n(k) . 
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B. The cluster expansions of the one-body density matrix 

The trial wave function 'J, which describes a correlated nuclear system, can be written as (e.g. |3fi|| ) 

* = T<S>, (8) 

where $ is a model wave function which is adequate to describe the uncorrelated A-particle nuclear system and T 
is the operator which introduces SRC. <E> is chosen to be a Slater determinant wave function, constructed by single- 
particle wave functions. Several restrictions can be made on the model operator T J36|,[3^]. In the present work T is 
taken to be of the Jastrow-type [|| 

A 

where /(r^ ) is the state-independent correlation function of the form 

f( nj ) = l-exp[-P(r l -r J ) 2 }. (10) 



1. Factor cluster expansion of Iwamoto-Yamada 



In the factor cluster expansion of Iwamoto-Yamada (FIY) the OBDM takes the form |T(| 

P fiy(t,t') = N[(O rr ,) 1 - 2 (r,r', gl ) - 2 (r,r',g 2 ) + 2 (r, r', g 3 )], (11) 
where N is the normalization factor, and the terms (O rr ')i and 2 (r, r', g;) (I — 1, 2, 3) have the general forms 

(O rr ,)i = PSD (r,r') = -Y^rintfl + l)<f>* nl (r)cl> n t(r')P l (cosL> rr ,), (12) 
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and 



where 



2 (r,r',g;) = J g ; (r,r',r 2 )[ j os_ D (r,r')psn(r 2 ,r 2 ) - psD{r, y 2 )psd{^2, r')]dr 2 , 



(13) 



gi(r,r',r 2 ) = exp[-/?(r 2 + r 2 2 )\ exp[2/3rr 2 ], g 2 (r,r',r 2 ) = g 1 (r',r,r 2 ), 

g 3 (r, r', r 2 ) = exp[-/3(r 2 + r' 2 )} exp[-2/3r 2 ] exp[2/3(r + r')r 2 ]. (14) 
The term 2 (r, r',gj), performing the spin-isospin summation and the angular integration, takes the general form 

2 (r,r', a )=4 ]T 

x 44; : ; : : : ; r ^r,r^g0-E(^0|fc0) 2 ^;; fc (r,r^g ^ ) . (ir,) 



where 



k=0 



<Kf (r,r',gi) = ^< (l (r) 4>n 3 i 3 (r') exp[-/?r 2 ] P h (co S u; rr ,) 



x / 0r* 2 ; 2 ( r 2)0r l4 ;4( r 2)exp[-/3r 2 ] i fc (2/3rr 2 )r 2 dr 2 , 



(16) 

Jo 

and the matrix element A™ 3 j 3 ™ 4 ; 4 ' fe (r, r', g 2 ) can be found from ( fL6| ) replacing r <-» r' and fi]7i <-* 713^3 while the 
matrix element corresponding to the factor g3 can be found from ( Jlq ) replacing the factors exp[— /3r 2 ], P; 3 (cosw rr ') 
and ife(2/3rr 2 ) by the factors exp[— (3(r 2 + r' 2 )], flf il3 (uj rr i) and ifc(2/3|r + r'|r 2 ) respectively Q. In the expressions 

of the matrix elements A™ 3 ^™ 4 ^' fe (r, r', g;), ik(z) is the modified spherical Bessel function and the factor f2^ 3 (o; rr ') 
depends on the directions of r and r'. 
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2. Factor cluster expansion of Aviles, Hartogh and Tolhoek 



In the factor cluster expansion of Aviles, Hartogh and Tolhoek (FAHT), truncated at the two-body terms, the 
OBDM takes the form (details of the calculations are given in Appendix I) , 



p F AHT(r,r') = -j(O rr >)i 



(A-l)(O r 



2 (r, r', gl ) - 2 (r, r', g 2 ) + 2 (r, r',g 3 ) 1 



A(A-l) - /[Oa(r,r,gi) + 2 (r,r,g 2 ) - 2 (r, r, g 3 )]dr A 



(17) 



where (O rr ')i and 2 (r,r',g;) are given again by Eqs. ( |12| ) an d respectively. The FAHT expansion has the 
advantage that the normalization is preserved term by term. 



3. Low order approximation 

In the low order approximation (LOA) of Gaudin et al Jl3| the Jastrow wave function ^ of the nucleus was expanded 
in terms of the functions g — f 2 (rij) — 1 and h — /(r»j) — 1 and was truncated up to the second order of h and the first 
order of g. This expansion contains one- and two-body terms and a part of the three-body term which was chosen so 
that the normalization of the wave function was preserved. In LOA the OBDM takes the form O,0,0| 



p LO A(r, r') = ^[(O rr ,)i - 2 (r, r', gl ) - 2 (r, r', g 2 ) + 2 (r, r', g 3 ) + 2<3 3 (r, r', 0) - 3 (r, r', 2/3)], (18) 

where (O rr <)i and 2 (r, r', g;) are given again by Eqs. ([lj) and ( |l3| ) respectively and the three-body term 3 (r, r', z) 
(z = f3, 2/3) has the form 

3 (r,r',z) = y*g(r 2 ,r 3 ,z)psD(r, r 2 )[psD(r2,r')/OsD(r3,r3) -/3SD(r 2 ,r 3 ) j osi)(i" 3 ,r / )] dr 2 dr 3 , (19) 

where 

g(r 2 , r 3 , z) = exp[-z(rj + r\ - 2r 2 r 3 )]. (20) 
The term 3 (r, r', z), performing the spin-isospin summation and the angular integration, takes the general form 

3 (T,r',z)=4 ^ VnikVn ] l j '>ln k l k (2k + 1) X 



TLili ,71 j Ij iTlfclfc 

4(2l k + l)6 hh Al^:^ {v,v',z) - {21, + l)5 hlk ^ (^OIjO | k'^A^^f (r, r', 

fc'=0 



where 

KfhZhZt k ( r > r ''^ = TZ^U^n^Phi^^,) x / C 2i2 (^)^ 5i5 (r 2 )exphzr 2 2 ]r 2 2 dr 



(21) 



47T 

/•oo 

x / C 3 ; 3 ( r 3) 0„ 6 i 6 (r 3 )exp[-zr^] exp[2zr 2 r 3 ] r^dr 3 . (22) 
Jo 

Expressions ([[2]), ( |l5| ) and (|2|) were derived for the closed shell nuclei with N = Z where r\ n i is or 1. For the 
open shell nuclei (with = Z) we use the same expressions where now < r\ n i < 1. The normalization is preserved 
for the closed shell nuclei in all the expansions. In the case of the open shell nuclei the normalization is preserved (in 
the above formalism) for FIY and FAHT expansions. In the case of LOA, in which the number of particles is also 
conserved |p^| , particular attention has to be paid in each open shell nucleus. 

It is noted that the general expressions of the two- and three-body terms of the density matrix given by Eqs. (|l5| ) 
and ( ^l| ) are also valid for the expansions of the DD, FF and MD. The only difference is the expressions of the matrix 
elements A which have to be used. For the DD they are found from ( |lq ) putting r' = r, while the ones of the FF 
follow from Eq. (||) replacing p(r) by A(r, r) and for the MD they follow from Eq. (Q) replacing p(r, r') by A(r,r'). 

In the case when the model wave function $ is constructed from harmonic oscillator (HO) wave functions, analytical 
expressions of the various terms of the DD, FF and MD for any N = Z s-p and s-d shell nuclei can be found for FIY 
and FAHT while in the case of LOA analytical expressions of the closed shell-nuclei in the same region can be found. 
These expressions which depend on the HO parameter b and the correlation parameter /3 are given in Refs. |l^ [Tt]] 
for FIY and LOA while the ones for FAHT can be found easily from the other expansions. 
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III. RESULTS AND DISCUSSION 



The three expansions, mentioned in Sec. II, have been used for the analytical calculations of the DD, MD and 
charge FF as well as for the calculation of the information entropy sum defined by Eq. (||). The HO parameter b 
and the SRC parameter (3 in the three cases have been determined, for each nucleus separately, by a least squares fit 
to the experimental charge FF as in Ref. |L5l (using the same expression for \ 2 ). The center-of-mass correction has 
been taken into account by a Tassie-Barker factor [p9| while those for the finite proton size and the Darwin-Foldy 



relativistic correction through the Chandra and Sauer approximation |40|. They are not taken into account in the 
calculations of DD and MD to obtain the information entropy sum (and in the plots of MD). 

The variation with A of the best fit values of the parameters b and (3 for each of the three expansions is shown in 
Fig. 1 where b and (3 versus the mass number A have been plotted for various s-p and s-d shell nuclei. It is seen that 
these parameters have the same behaviour in FIY and FAHT expansions. In the case of LOA expansion, which has 
been used only for 4 He, 16 and 40 Ca the variation of the parameters seems to be the same. From Fig. lb it is seen 
also that the SRC parameter (3 has larger values in the open shell nuclei ( 12 C, 24 Mg, 28 Si and 32 S) than in the closed 
shell ones, indicating that there should be a shell effect in the case of closed shell nuclei. 

In this work we compare different expansions on the example of MD for closed and open shell nuclei. The reason 
for this is that the high- momentum component of n(k) is very sensitive to the extent to which nucleon correlations 
are accounted for in a given correlation method and in various approximations. The effect of different expansions on 
the form factors can be seen comparing the values of x 2 f° r the various expansions and nuclei. 

The MD for the closed shell nuclei 4 He, 16 and 40 Ca, calculated with the best fit values of the parameters 
and for the three expansions, are shown in Fig. 2. It is seen that the inclusion of SRC increases considerably the 
high momentum component of n{k). It has the same slope up to 2 fin - for the three expansions. In the region 
2 fm _1 < k < 5 im~ the slope seems to be a little different. FIY gives a larger contribution in the high momentum 
component than FAHT and LOA which give the same contribution in this region. The same behaviour of n(k) has 
been observed in the open shell nuclei as can be seen from Fig. 3. Here we would like to note that in general, a more 
realistic description of MD requires the inclusion of tensor correlations in the theoretical scheme. 

In the previous analysis, the nuclei 24 Mg, 28 Si and 32 S were treated as Id shell nuclei, that is, the occupation 
probability of the 2s state was taken to be zero. The formalism of the expansions FIY and FAHT has the advantage 
that the occupation probabilities of the various states can be treated as free parameters in the fitting procedure of 
the charge FF. Thus, the analysis can be made with more free parameters. For that reason we considered, as in Ref. 
|fl6f the cases FIY* and FAHT* in which the occupation probability r\2s of the nuclei 24 Mg, 28 Si and 32 S was taken to 
be a free parameter together with the parameters b and (3. We found that in both expansions the x 2 values become 
smaller, compared to those of cases FIY and FAHT and the A dependence of the parameter (3, as can be seen from 
Fig. lb, is not so strong as before. Also the values of rj2 S found in the fit and the values of rjid found through the 
relation rji^ = [(Z — 8) — 2rj2 S ]/10, are very close for both expansions in each nucleus. 

Our best fit values of the parameters and the values of \ 2 for the various nuclei under consideration and for the 
three expansions as well as for the HO case (that is when SRC are not included) are shown in Table I. From the 
values of x 2 we conclude that the three expansions give similar values of x 2 - The FIY and FAHT expansions have 
almost the same x 2 values. They differ less than 2% in the two expansions in each nucleus. In most cases the x 2 
values corresponding to FIY (or FIY*) are smaller. There are two cases ( 12 C and 28 Si) when the FAHT or FAHT* 
expansion gives smaller x 2 value and one case ( 16 0) when LOA gives smaller x 2 value. 

In addition, we verify the information-theoretic criterion for comparing the quality of the three expansions. It is 
seen in Table I that almost in all cases, the larger the S the smaller the \ 2 . Both methods of comparison (S and x 2 ) 
show that the FIY (or FIY*) expansion is better than the FAHT and LOA for 4 He, 24 Mg, 32 S and 40 Ca. For 16 the 
LOA is the best. There are only two exceptions to this rule i.e. in 12 C for cases FIY and FAHT and in 28 Si for cases 
FIY* and FAHT*. In 12 C x 2 is smaller in FAHT and we expect S to be larger than in FIY while in 28 Si x 2 is smaller 
in FIY* and we expect S to be larger than in FAHT*. These are two exceptions to our rule. It should be noted also 
that in these two exceptions the difference in the x 2 values for the two expansions in both nuclei is less than 1%. 

Finally, we consider the so-called "healing" or "wound" integrals, denoted here as w 2 nl [ j37pl] ] for the various states 
of the relative two-nucleon motion, pertinent to the closed shell nuclei of Table I and in each case, that is in each of the 
cluster expansions FIY, FAHT and LOA. The values of these integrals express in a way the "amount of correlations" 
introduced to each state of the relative two-nucleon motion. The healing integrals (for a state independent correlation 
function /(r), such as the one given by (0)) are defined as follows 



■■>rd ~ I \tpni(r) ~ (t>ni(r)\ dr, (23) 
Jo 

where <j) n i(r) is the (normalized to unity), uncorrelated (HO) radial relative wave function and ip n i(r) the correspond- 
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ing, normalized to unity, correlated one: ip n i(r) — N n if(r)<p n i(r), where N n i, the normalization factor of ipni^), is 
given by 

-1/2 

(24) 

It is interesting to note that with the correlation function (|T(]) the healing integrals can be calculated analytically 
for every state nl. Some details are given in Appendix II. As one expects, these integrals depend on both, the HO 
parameter b and the correlation parameter [3. We may note, however, that their dependence on them is only through 
the dimensionless product y — 2f3b 2 (see expression ( |44| ) of Appendix II). 

In Table II the values of the parameters b, (3 and y = j3b 2 for each closed shell nucleus and cluster expansion 
considered, are displayed along with the corresponding values of for certain relative state in the s-p and s-d closed 
shell nuclei. It is seen from the results in this table that the values of u>^ ; , for each of the relative states (nl) involved 
in each nucleus, are smaller when is obtained with the FIY expansion and larger when obtained with the LOA. 
Furthermore, for each nucleus and expansion the values of w^i of the nodless (n=0) states decrease as the value of 
I increases, the correlations having less effect to these higher /-states, because of the existing centrifugal (repulsive) 
term of the HO potential. The values of w^q increase when n — 1 or n = 2 in comparison with those of Wq . 

IV. SUMMARY 

In the present work, a systematic study of the effect of SRC on one-body properties of sp and sd shell nuclei has 
been made evaluating three different cluster expansions. The HO parameter b and the SRC parameter (3 have been 
determined by a least-squares fit to the experimental charge FF. 

The comparison of the three expansions on the example of the MD and the FF shows that they can be considered 
as equivalent expansions. It is found that, when the calculations are made with the best fit values of the parameters, 
these expansions reproduce the diffraction minima of the FF in the correct place and they give similar MD for all the 
nuclei we have considered. The inclusion of SRC increases considerably the high momentum component of n(k). 

The FIY and FAHT expansions have been used both for closed and for open shell nuclei while the occupation 
probabilities can be treated as free parameters together with the parameters b and f3 in the fitting procedure of the 
FF. In LOA such calculations are in progress. 

In addition, the information entropy sum has been calculated according to the three methods compared in the 
present work. It was found almost in all of the numerous cases (different expansions and nuclei), that the larger the 
S, the smaller the x 2 ■ That is S could be used as a criterion for the quality of a given nuclear model. We found only 
two exceptions to this rule. In these two exceptions the difference of the % 2 values is less then 1%. 

Finally, attention was paid to the "healing" or "wound" integrals w 2 nl of the relative two nucleon states. A convenient 
analytic expression of with correlation function (|l(]) was derived for any relative state nl. Their values were 
computed in a number of states with that expression and were also discussed. 

V. APPENDIX I 

In this appendix, we give some details about the FAHT expansion. We define the correlated wave function as 

A 

* = n/(r l ,r J )$, (25) 

i<j 

where f(ti, Vj) is the Jastrow correlation function and $ is a Slater determinant wave function. To built up the cluster 
expansion, we start, following Ref. p0|] , from the A-body integrals Ja(A) defined as 

A A A 

^(A)^ ]t (4...fein/( r ^ 3 )°i(v o!(A) n/( r ^;)i4-i)« > ^ 

h—tA i<j i<j 

where the sum over the states ii, ia> • • • >M has no restrictions and extends over all one-particle states and a stands 
for the antisymmetrization. The operators Oi(^4) and 02(^4) have the forms 



N n l = 



f(r)4> 2 nl (r)dr 
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1 (A) = l[S(r i -r' i ), 

lli=i °V r » r J i=l 



The OBDM pfaht(t,t'), normalized to A, is denned as 

PFAHT{r,r') = 

We introduce the n-body integrals J n (A) defined as 

1 



dlnJ A (A) 



dA 



A=0 



J„(A) 



A(i4-l)---(A-n + l) 



E 



n/(r„r i )0 1 (n)e AO! 'Wn/(r' i y j )|^ 



ii • ■ ■ Vi n /a 



H...l n i<j i<j 

The cluster integrals 3„ (n — 1, 2, . . . A) are defined through the successive application of the equation 



Jn = U 



\5- 



5 n , 7i — 1,2,..., A. 



k=l 



For example, for n = 1 and n = 2 it gives 



3?i = Ji, 3 2 = 



/2 



The last of Eqs. (p0[), which corresponds to n = A is the quantity we are interested in 



J A =U 



,(0 =9 (f) 9 UU0 



If the factor-cluster expansion is limited to the two-body term (assuming that the remaining cluster inte; 
equal to unity j2(J]), then 

T ^, (^( 1 )c^( 2 ) 
J A « «i 3 2 ■ 



From Eqs. (|28j) and (|33|) we have 
/) M if T (r,r') = 





' 1 


dJi" 






' 1 


dJ 2 


(0 






+ 

A=0 


(0 


.^2 


dA 



- 2- 



A=0 



where 



and 



il=l 



(27) 

(28) 
(29) 

(30) 

(31) 

(32) 
;rals are 

(33) 
(34) 

(35) 
(36) 



After some algebra we obtain 
Ji(0) = 1, 



Ja(0) 
dJi 



1 



A(A - 1) 



A{A-1)- / [0 2 (r,r,. gi ) + 2 (r,r,.g 2 )-0 2 (r,r, 5 3)]dr 



dA 
dJ 2 
dA 



A=0 



1 

T 



(Or 



x=o A{A-l) 



{A — l)(O rT ./)i - 2 (r,r', 5l ) - 2 (r,r', 52 ) + 2 (r,r', ff3 ) 



(37) 
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where the terms (O rr >)i and 2 (r, r',gi) have been defined in Sec. II. 
Finally, the Pfaht{ y i r ')i normalized to unity, becomes 



Pfaht i"') = -^(O 



(38) 



-(A-l) 



(A-lXOg^ - Q 2 (r,r', gl ) - Q 2 (r,r',g 2 ) + Oa^rVga) 1 
A(A - 1) - /[0 2 (r, r, 5l ) + 2 (r, r, 52 ) - 2 (r, r, 53 )]dr ^ 



O, 



VI. APPENDIX II 



The healing integral defined by(p3|) is written as follows fl4l| 



where 



u 2 nl = 2 [1 + N nl (I nl (b,P)-l)}, 



I nl (b,P)= / exp[-/3r 2 ]^(r)dr, 
Jo 



(39) 



(40) 



and the normalization factor N n i is given by (g4|) . This factor can be easily expressed in terms of the integrals I n i (b, /3) 
and I n i(b, 2/3) by means of expression (UG) 



N nl - [1 - 2I nl (b,f3)+I nl (b,2{3)]- 1/2 . 



(41) 



Thus, the analytical calculation of any healing integral is reduced to the calculation of two integrals of type 
(|40|). The expression of I n i{b, 2/3) follows immediately from the expression of I n i(b, (3). 
We use the general expression of the radial HO wave function (normalized to one as J °° (j^idr = 1) in the form 



(f>ni(r) 



2n\ 



T(n + l + ±)b r 



1 -«+* fr 2 



— ) L n 2 ( p ) exp 



2bl 



(42) 



where b r is the HO parameter of the relative motion, which is related to the usual HO parameter b by b r = \[2b 
{b= {h/muj) 1 / 2 ). 

Substituting expression @ into @ and using the transformation r 2 jb 2 — £, /„; is written 



Inl(b,P) 



V[n + l + 3/2] J 



-(i+y)?^+i/2 



(43) 



where y = fib 2 r = 2f3b 2 . 

Using formula 13 of §7.414 of Ref. Gl I nl takes the form 



l nl (b,0) = (y~ir(y + l)— l -^P^o) (V±^ 



(44) 



where Pi° 1,a2 '(z) the Jacobi polymomials. These may be easily expressed in terms of the Hypergeometric function 



(see e.g. §8.962 of Ref. [|42|). In the case of the nodless states (because P^ ai ' a2 \z) = 1) /„; takes the simple form 



I 0l (b,P) = (y + l)- l -V 2 , 



(45) 



By substituting (3 — > 2/3, the expression of I n i(b, 2/3) follows immediately and therefore the analytic expression of 
the w^i by means of the formulae ( |39| ) and (|J). It is thus clear that the healing integral w^i for any state depends 
on correlation parameter /3 and the HO one, only through the product y — 2(3b 2 . The expressions of for the lower 
n— states follow also very easily. 
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FIG. 1. The harmonic oscillator parameter b (a) and correlation parameter f3 (b) versus the mass number A for the expansions 
FIY, FAHT and LOA. Cases FIY* and FAHT* (open circles and squares respectively) correspond to the case when the 
occupation probability 7/2s is treated as a free parameter. 




k [fm" 1 ] k [fm" 1 ] k [fm" 1 ] 



FIG. 2. The momentum distribution of the closed shell nuclei in the three expansions, FIY (solid line), FAHT (dash line) 
and LOA (dot line). The normalization is J n(k)dk = 1. 
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k[fm' 1 ] k [fm 1 ] 



FIG. 3. The momentum distribution of the open shell nuclei in the case of the FIY expansion (solid line) and the FAHT 
expansion (dash line). The normalization is as in Fig. 2. 
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TABLE I. The values of the parameters 6 (in fm) and j3 (in fm~ 2 ), the \ 2 > the RMS charge radii {r 2 ^ 1 ^ 2 (in fm), of the 
mean kinetic energy per nucleon (T) (in MeV) and the nuclear information entropy in position- (5V) and momentum-space 
(Sk) and the sum of them S for various s-p and s-d shell nuclei. The various cases have been ordered according to increasing 
values of \ 2 ■ For the various cases see text. 



Nucleus 


Case 


6 


f3 


x 2 


{r 2 ch ) 1/2 


(T) 


Sr 


s k 


S 


4 He 


FIY 


1.1732 


2.3127 


3.50 


1.623 


29.904 


9.978 


5.985 


15.963 




FAHT 


1.1661 


1.9092 


3.70 


1.621 


29.048 


9.943 


6.013 


15.955 




LOA 


1.1605 


1.6584 


3.88 


1.620 


28.543 


9.917 


6.034 


15.951 




HO 


1.4320 


oo 


30.94 


1.765 


15.166 


11.632 


3.014 


14.646 


12 C 


FAHT 


1.5204 


2.4683 


90.19 


2.427 


24.779 


31.455 


1.989 


33.444 




FIY 


1.5190 


2.7468 


90.87 


2.426 


25.580 


31.436 


2.142 


33.578 




HO 


1.6251 


oo 


176.54 


2.490 


17.010 


32.714 


-2.2484 


30.465 


16Q 


LOA 


1.6387 


1.8825 


115.50 


2.674 


23.006 


42.083 


-4.393 


37.690 




FIY 


1.6507 


2.4747 


120.19 


2.680 


23.614 


42.237 


-4.557 


37.680 




FAHT 


1.6554 


2.2097 


122.49 


2.684 


22.518 


42.313 


-4.939 


37.374 




HO 


1.7610 


oo 


199.45 


2.738 


15.044 


43.655 


-10.667 


32.988 


OA-k .1- 

24 Mg 


FIY* 


1.7473 


2.4992 


140.37 


3.064 


24.614 


63.532 


-14.334 


49.198 




FAHT* 


1.7468 


2.1833 


140.40 


3.064 


23.742 


63.536 


-14.603 


48.933 




FIY 


1.8103 


4.2275 


177.51 


3.095 


21.109 


64.452 


-19.228 


45.224 




FAHT 


1.8120 


4.1322 


177.91 


3.096 


20.818 


64.483 


-19.410 


45.073 




HO 


1.8496 


oo 


188.01 


3.117 


16.162 


65.124 


-23.429 


41.695 


28 Si 


FAHT* 


1.7773 


2.1193 


103.39 


3.184 


24.184 


72.901 


-20.844 


52.057 




FIY* 


1.7774 


2.4440 


103.47 


3.184 


25.205 


72.888 


-20.438 


52.450 




FIY 


1.8236 


3.0020 


126.33 


3.216 


22.933 


73.889 


-24.115 


49.774 




FAHT 


1.8279 


2.8372 


127.84 


3.219 


22.110 


73.987 


-24.645 


49.342 




HO 


1.8941 


oo 


148.28 


3.257 


16.099 


75.288 


-32.022 


43.266 


32q 


FIY* 


1.8121 




1 1SIS 1 1 


3.282 


94 Q1 fi 


89 i nn 

O^. -LOO 


98 "H4^ 


OO. ( oo 




FAHT* 


1.8131 


2.3358 


166.31 


3.283 


23.961 


82.129 


-28.827 


53.302 




FIY 


1.9368 


3.0659 


304.96 


3.443 


20.867 


86.921 


-36.707 


50.214 




FAHT 


1.9417 


2.9585 


306.46 


3.446 


20.252 


87.045 


-37.316 


49.729 




HO 


2.0016 


oo 


320.45 


3.483 


14.878 


88.361 


-44.881 


43.480 


40 Ca 


FIY 


1.8660 


2.1127 


160.44 


3.516 


26.617 


101.501 


-42.710 


58.791 




FAHT 


1.8685 


1.7397 


161.13 


3.517 


24.643 


101.558 


-44.172 


57.387 




LOA 


1.8164 


1.7404 


188.36 


3.397 


25.586 


97.611 


-42.121 


55.490 




HO 


1.9453 


oo 


229.32 


3.467 


16.437 


100.987 


-58.709 


42.278 
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TABLE II. The values of the parameters b (in fm), fi (in fm 2 ) and y = /3b 2 and the values of the healing integral w^i for 
various states and for the closed shell nuclei 4 He, le O and 40 Ca and the three expansions FIY, FAHT and LOA. 



Nucleus 


Case 


b 


P 


y = f3b z 


2 


2 


2 




2 

W>10 


2 

">03 


4 He 


FIY 


1.1732 


2.3127 


3.1832 


0.01874 














FAHT 


1.1661 


1.9092 


2.5961 


0.02450 














LOA 


1.1605 


1.6584 


2.2335 


0.02971 












16 o 


FIY 


1.6507 


2.4747 


6.7431 


0.00664 


0.00024 


8.6 10" 


-6 


0.00925 






FAHT 


1.6554 


2.2097 


6.0554 


0.00773 


0.00031 


1.2 10~ 


-5 


0.01069 






LOA 


1.6387 


1.8825 


5.0552 


0.00996 


0.00048 


2.3 10~ 


-5 


0.01359 




40 Ca 


FIY 


1.8660 


2.1127 


7.3563 


0.00586 


0.00020 


6.4 10" 


-6 


0.00821 


2.1 10~ 7 




FAHT 


1.8685 


1.7397 


6.0738 


0.00770 


0.00031 


1.2 10" 


-5 


0.01065 


4.8 10~ 7 




LOA 


1.8164 


1.7404 


5.7421 


0.00833 


0.00035 


1.5 10" 


-5 


0.01148 


6.2 10~ 7 
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